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^0 ■ Aims. In this paper we aim to determine the structure on 100 AU scales of the massive young stellar object W33A, using interferometric 
pH ' observations in the mid-infrared. This emission could be due to a variety of elements, for example the inner protostellar envelope, outflow 
f») cavity walls, dusty or gaseous accretion disk. 
I Methods. We use the Unit Telescopes of the VLT Interferometer in conjunction with the MIDI instrument to obtain spectrally dispersed 
O ' visibilities in the Af-band on 4 baselines with an angular resolution between 25 and 60 milli-arcseconds (equivalent to 95 and 228 AU at 3.8 kpc). 
The visibility spectra and spectral energy distribution are compared to 2D-axi-symmetric dust radiative transfer models with a geometry 
including a rotationally flattened envelope and outflow cavities. We assume an 07.5 ZAMS star as the central source, consistent with the 
observed bolometric luminosity. The observations are compared to models with and without (dusty and gaseous) accretion disks. 
Results. The visibilities are between 5% and 15%, and the non-spherically symmetric emitting structure has a typical size of 100 AU. A 
satisfactory model is constructed which reproduces the visibility spectra for each (u,v) point. It fits the A^-band flux spectrum, the mid-infrared 
slope, the far-infrared peak, and the (sub)mm regime of the SED. It produces a 350 jim morphology consistent with the observations. 
Conclusions. The mid-infrared emission of W33A on 100 AU scales is dominated by the irradiated walls of the cavity sculpted by the outflow. 
The protostellar envelope has an equivalent mass infall rate of 7.5 10~ 4 M Q yr _1 , and an outflow opening angle of 29 = 20°. The visibilities 
rule out the presence of dust disks with total (gas and dust) masses more than 0.01 M Q . Within the model, this implies a disk M acc of less than 
1.1 10~ 7 (or/0.01) M Q yr -1 , where a is the viscosity of the Shakura-Sunyaev prescription. However, optically thick accretion disks, interior to 
the dust sublimation radius, are allowed to accrete at rates equalling the envelope's mass infall rate (up to 10 3 M Q yr ) without substantially 
■ affecting the visibilities due to the extinction by the extremely massive envelope of W33A. 
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5_i 1. Introduction growing number of direct and indirect pieces of evidence point 



to the presence of large-scale (10,000 AU) rotating toroids (e.g. 

Circumstellar accretion disks are an essential element in the Bekran et al 2004i 2005; Fumya et al 2008)> but also t0 in _ 

formation of stars. The relatively long time scales involved in ^g^y disk-like structures on scales of several 100 AU in 

low-mass star formation (SF) and the proximity of the regions a smal] number of comparatively nearby examp i e s (e.g. Hoare 

where they form allow the disks to be imaged at millimeter 2006; j im6nez -Serra et al. 2007; for a review see Cesaroni et 

(Dutreyetal. 1996) , near-infrared (near-IR; Padgett et al. 1999) a j 2001) 
or even optical wavelengths (Grady et al. 1999) and studied in 

great detail. The short time scales on which high-mass stars are The reality of a disk accretion scenario in massive star for- 

formed and the distance of massive SF regions amongst others mation can be assessed by observations of massive stars during 

render the study of the accretion process particularly difficult, the main mass gathering phase at significant angular resolution. 

Whether or not a high-mass star grows in mass through the ac- Prime examples of such young massive stars constitutes the 

cretion of matter delivered to the stellar surface by means of class of object referred to as massive YSO (MYSO), but also as 

a circumstellar disk is therefore still left unanswered. Recent high-mass protostellar object, and BN-object. The class is char- 

3-D radiation hydrodynamic simulations demonstrate that ac- acterised by luminous (> 10 4 L ) infrared objects with spectral 

cretion through a disk can continue in the presence of strong energy distributions that peak around 100/vm. The luminosity 

radiation pressure (Krumholz et al. 2009). Observationally a indicates a single ZAMS star that would be able to ionise its 

surroundings, however too little or no radio emission is ob- 

Send offprint requests to: W.J. de Wit, e-mail: served. This absence has been used as an argument to claim 

w.j.m.dewit@leeds.ac.uk on-going mass infall onto the central object. As a result the 

* Based on observations with the VLTI, proposal 381.C-0602 star is extended and thus relatively cool (Hosokawa & Omukai 
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2009), and therefore unable to ionise its surroundings (Hoare & 
Franco 2007). Alternatively, high mass infall rates (Walmsley 
1995) or gravitational entrapment (Keto 2002) could quench 
the development of an Hn that dominates the radio contin- 
uum emission. These alternative scenarios invoke however very 
high emission measure (EM) Hn regions to make them opti- 
cally thick in the radio, e.g. Ti cm = 3 10~ 10 x EM(cirT 6 pc) 
(Osterbrock 1989). However, it is unlikely in this picture that 
the region would still be optically thick in the near-IR H i re- 
combination lines Bra- and Bry as these have line centre op- 
tical depths that are 10 5 ~ 6 times smaller than ri cm (Hummer 
and Storey 1987). This would imply that these objects should 
have strong near-IR Hi lines and yet they only show weak, 
broad lines indicating a stellar wind origin (Bunn et al. 1995). 
Nonetheless, the MYSO class has the potential to test the disk 
accretion scenario. The actual accretion region is probably 
found on AU to tens of AU scales near the central stellar object. 
This region is however buried by hundreds of A v of dust, which 
makes up the extended (~0. 1 pc) protostellar envelope (van der 
Tak et al. 2000). Only long-wavelength radiation is capable of 
carrying away information on the physical processes at play, 
limiting the attainable spatial resolution with single dish aper- 
tures. Interferometers can overcome the problem of limited res- 
olution and the near and mid-infrared wavelength regions cur- 
rently deliver the highest angular resolution by means of e.g. 
the Very Large Telescope Interferometer (VLTI). 

The inner parts of the dusty MYSO protostellar enve- 
lope is heated to a few 100K and emits strongly at mid-IR 
wavelengths. The 10 /mi emission is usually unresolved at the 
0.3" resolution of single 8 meter dish telescopes, but it is 
partially resolved at 24.5/mi (de Wit et al. 2009). Here, we 
present new observations made of the MYSO W33A with the 
MID-infrared Interferometric (MIDI) instrument at the VLTI. 
MIDI is a two beam combiner that operates in the TV-band (8- 
13/an) and delivers spectrally dispersed interferometric observ- 
ables (see Leinert et al. 2003) . W33A is one of the most in- 
tensely studied MYSO thanks to its IR spectrum rich in ice fea- 
tures (Gibb et al. 2000). Its physical nature has received much 
less attention. The object has an IRAS luminosity of 10 5 L G 
(Faundez et al. 2004) for a kinematical distance of 3.8 kpc 
(Bronfman et al. 1996). The object shows weak, compact and 
optically thick radio continuum emission (Rengarajan & Ho 
1996; van der Tak & Menten 2005) and broad (-100 kms" 1 ), 
single-peaked Hi recombination emission consistent with an 
ionised stellar wind origin (Bunn et al. 1995). Modelling of 
the protostellar envelope with spherical dust radiative transfer 
(RT) model at constant density by Giirtler et al. (1991) required 
the presence of an inner dust-free cavity of 135 AU (35milli- 
arcsecond at 3.8 kpc), a size that would be well resolved by 
MIDI. In de Wit et al. (2007, hereafter paper I) we presented 
a single MIDI baseline for W33A in order to probe this hun- 
dred AU scale region. We found using 1-D modelling an inner 
dust-free radius at 7 milli-arcsecond (26.5 AU), corresponding 
to the dust sublimation radius. We concluded that the dispersed 
MIDI visibilities and SED can be reproduced by a spherical 
dust model provided that the density law is relatively flat. In 
this paper we present three additional VLTI/MIDI observations 
of W33A with baseline position angles close to perpendicular 



Table 1. MIDI observations of W33A executed in high-sens 
mode at a spectral resolution of R = 30. Seeing is the one mea- 
sured by the DIMM. Configuration D was already presented in 
de Wit et al. (2007). 
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2008-06-21 
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2008-04-22 
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1.08 
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2005-09-16 


U2-U3 


45.58 


47.31 
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and with different projected lengths. This allows us to constrain 
the distribution of the emitting material further using 2-D axi- 
symmetric RT models. 

We describe and discuss the new and archived MIDI obser- 
vations in Sects. |2] and |3] The basic features of the RT model, 
our modelling approach and the final preferred model are ex- 
plained in Sect. [4] The implications for W33A and their effect 
for MYSOs in general are presented in Sect. [5] We conclude 
our work in Sect.|6] 

2. Observations and data reduction 

2.1. MIDI observations 

W33A was observed with the VLTI and the MIDI instrument 
in service mode on four different occasions using the UTs as 
aperture stations (see TableQ}. We aimed at covering differ- 
ent position angles on the shortest possible UT baselines. The 
baselines are presented in Fig.Q]projected onto the near-IR re- 
flection nebula of W33A as imaged by UKIDSS (Casali et 
al. 2007; Warren et al 2007; Lucas et al. 2008). The shortest 
VLTI baselines are attained with the ATs but W33A cannot 
be observed in this mode due to its relatively low brightness. 
The obtained angular resolution for the data set presented here 
stretches from 25 to 60 milli-arcseconds. Detailed descriptions 
of the MIDI observation procedure are given in Przygodda et 
al. (2003), Chesneau et al. (2005) and Leinert et al. (2004). 
The W33A observations were executed in the so-called High- 
Sens MIDI mode, which uses all the incoming light for beam 
combination and fringe tracking. The incoming beams are in- 
terfered producing two complementary interferometric chan- 
nels that have by definition a phase difference of n radians. A 
measurement of the uncorrected flux spectrum to obtain final 
visibilities is done immediately after the interferometric obser- 
vation. The high-sens mode is advantageous when observing 
faint targets or small visibilities, however the accuracy of the 
final visibility spectrum is limited by the sky brightness vari- 
ation between the interferometric and photometric measure- 
ment. This can amount to 10-15% uncertainty, depending on 
the atmospheric stability. During the nights, the unresolved in- 
terferometric calibrator star HD 169916 was observed in or- 
der to determine the instrumental visibility and correct for it. 
We also inspected all other observed calibrators from differ- 
ent programs for each night (a total of 18 objects), which al- 
lowed us to estimate the stability and accuracy of the instru- 
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mental visibility. We found it to be uncertain by not more than 
5%, much smaller than the uncertainty in the determined pho- 
tometric measurement (20%, see next section). The calibrator 
star provides also an approximate flux calibration (Cohen et al. 
1999) in addition to the visibility calibration. A prism with a 
spectral resolution of 30 was employed in order to disperse the 
incoming beams. 

The correlated flux spectra are estimated by the MIA+EWS 
software package (version 1.6; see Jaffe 2004; Kohler 2005). 
Interferograms are coherently added in order to maximize the 
signal to noise. First, fringe spectra have to be corrected, be- 
cause they are affected by the instantaneous atmospheric and 
instrumental piston. The fringe scan is used to estimate the 
group delay due to the atmosphere. Removing the atmospheric 
and the (known) instrumental group delays constitutes a cor- 
rection to the dispersed fringe signal, and straightens the dis- 
persed fringe spectra, i.e. the phase is independent of wave- 
length. Next, the phase offset due to varying water refraction 
index between the time of recording of the fringe spectra has to 
be accounted for. In principle all spectra can then be added to 
a final fringe spectrum and the amplitude of the power at each 
frequency estimated, which will give the correlated flux. 

The total flux spectra are taken in a chopping mode and 
each observation results in two spectra as MIDI splits the 
beams. The spectra are extracted from the raw data in the EWS 
mode. An estimate of the sky background is made and sub- 
tracted. The area corresponding to the mask with which the cor- 
related spectra are extracted is used here. The counts are then 
simply summed in the y-direction. The final spectrum is the 
sum of the two geometric means of the four individual spectra: 
VA1 x Bl + VA2 x B2. This quantity is also what is obtained 
for the correlated flux after beam combination and thus ensures 
consistency in deriving visibilities. 

2.2. JCMT/HARP-B observations 

W33A was observed on the night of 5 June 2007 in 12 CO, 
13 CO and C ls O using the HARP-B instrument on the James 
Clerk Maxwell Telescope. The 12 CO data have a resolution 
of 488 kHz, and the other two lines were rebinned to the 
same resolution, giving a velocity resolution of approximately 
0.4 km s -1 . The maps had a linear baseline removed. The core 
emission was presumed to be traced by the C ls O emission 
(which has line profile very close to a Gaussian). We used this 
as a model to trace the emission in the 12 CO and 13 CO maps 
that arises from the actual outflow. Full details of this work, 
which is part of a larger project studying outflows from the 
Red MSX Source Survey (see e.g. Urquhart et al. 2008), will 
be published separately in due course. 

3. Description of observational data of W33A 

3.1. N -band from MIDI 

The measured MIDI observables (flux, correlated flux) and de- 
rived visibilities of W33A are presented in Fig. [3] Each column 
represents a VLTI configuration and we have included the ob- 
servations already presented in paper I (configuration D). The 




Fig.l. UKIDSS K-band observations of W33A. Projected 
VLTI/MIDI baselines are indicated by white bars. 

top panels show the flux calibrated A^-band spectra. They show 
a variation of around 20% as expected for data taken in the 
observational high-sens mode. The errorbars represent two un- 
certainties that are added in quadrature: (1) the rms estimates 
of the flux and (2) the systematic error determined from a total 
of five W33A MIDI spectra (more than the standard one flux 
spectrum per interferometric measurement were taken.). The 
distinct feature of W33A in the A^-band is the extremely deep 
silicate absorption, the actual depth is not detected at the cen- 
tral wavelength where it drops to less than 1 % times the level 
of the pseudo-continuum. 

The second row of panels shows the correlated flux spectra. 
The dotted lines indicate the correlated noise level. It was esti- 
mated by processing that part of the detector directly adjacent 
to the W33A fringes, indicative of the background noise. The 
noise produces a signal a few times 0.1 Jy around 10 fim, de- 
spite the absence of incoming flux at the base of the silicate ab- 
sorption. Similar noise levels have been found previously (e.g. 
Jaffe et al. 2004; Matsuura et al. 2006). The errorbars of the 
W33A correlated flux spectra reflect the uncertainties due to 
spurious correlated noise level. We discard any signal below 
this noise level in further analysis. 

The third row of figures presents the resulting visibilities 
for the wavelengths not affected by low flux levels. The error- 
bars are the uncertainties from correlated and total flux added 
in quadrature. This entails that the systematic uncertainty inher- 
ent to the high-sens mode are taken into account in the visibility 
uncertainties. We find that W33A has visibility values between 
5% and 15% on the used baselines. The spectral shape of the 
visibilities reveal a clear trend of decreasing visibilities on the 
blue wing of the silicate absorption feature. The smaller vis- 
ibilities imply a larger emitting region, and therefore that the 
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Fig. 3. MIDI observations of W33A on four different occasions (see TableHJ. Top: Flux spectra. Middle: Correlated flux spectra 
(points with errorbars), and the detector noise level (dashed curve). Bottom: Derived visibilities. 



increasing optical depth due to the silicate absorption is due 
to structures on larger size scales than the pseudo-continuum. 
The visibilities corresponding to the two similar configurations 
(B and C) are practically the same, as are their correlated flux 
levels. It gives an independent indication of the stability of the 
set-up, and the relatively good weather. Configuration A shows 
a correlated flux level twice higher than the ones of configura- 
tion B and C, despite the relatively little difference in total flux. 
Given that the weather during the observations was stable, we 
believe that the difference in visibilities between the various 
configurations is real. If we translate the visibilities at 8 fim to 
Gaussian emission regions then the FWHM decreases from 
1 15 AU to about 95 AU rotating in P.A. from 15°(configuration 
A) to ~ 120° (B and C) and changing the baseline length ac- 
cordingly. 



3.2. The molecular outflow and near-IR reflection 
nebula 

The star formation activity in W33A produces an archetypical 
bipolar molecular outflow (Fig.|2]l. The south-eastern lobe con- 
stitutes the approaching, blue-shifted flow. Detailed millime- 
ter line maps which would provide constraints on the outflow 
opening angle and the inclination of the system are currently 
unavailable. However, outflow activity produces reflection neb- 
ulae particularly well traced by near-IR scattered light (see e.g. 
Tamura et al. 1991). The near-IR image of W33A already pre- 
sented in Fig. [T]is consistent with the CO map in position angle 
and orientation. We use this image therefore in order to con- 
strain these system parameters, which will be of further help 
in modelling the MIDI visibilities and SED later on. The full 
extent of the ~ 50" long /T-band nebula is shown in Fig.Q] The 
nebula has a position angle of ~ 145°, however on scales of 
order 1" patchy extinction distorts the morphology. A fainter 
extension can also be seen in the north-western direction. This 
lends support for the 145° position angle to be the actual orien- 
tation of the system. High-angular resolution IFU images pre- 
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Arc Seconds 

Fig. 2. 12 CO (3 - 2) map taken with JCMT/HARP illustrating 
the molecular outflow of W33A. The P.A. and orientation of 
the blue and red-shifted lobes of the outflow are consistent with 
the A'-band reflection nebula of Fig. Q] Coordinates are in J2000 
and they correspond to the UKIDSS near-IR source, and are in 
accord with the Capps et al. (1978) main IR source. 

sen ted in Davies et al. (2010) confirms the reflected nature of 
the ZT-band emission at the base of the outflow and position 
angle. The outer regions of the nebula may have contributions 
from shocked H2 emission as traced by extended 4.5 fim emis- 
sion in Spitzer IRAC images (Cyganowski et al. 2008). 

In Davies et al. (2010), it is also shown using spectro- 
astrometry of the Bry transition that the near-IR source is at 
the base of the bipolar outflow on scales of 200yU-arcseconds. 
The near-IR source corresponds most likely to the central ob- 
jects in W33A, indeed as was assumed during the execution 
of the MIDI observations. This source shows a Zf-band pro- 
file in the UKIDSS image that is broader than the unresolved 
stars in the nearby field. In the //-band the scattering nebula 
is much weaker (see Fig.|4|. Given that photons are more effi- 
ciently scattered at shorter wavelength the relatively dim nature 
of the //-band nebula argues for strong extinction in the line- 
of-sight towards W33A. Moreover, the central source is absent 
in the //-band. 

Most of the emission seen at Zf-band is dust scattered; con- 
tribution to the emission by molecular hydrogen gas is small as 
evidenced by intermediate resolution spectra, obtained as part 
of the RMS survejfl. In the //-band a foreground star is visible 
projected onto the nebula, located 2.5" south and 1.0" east of 
the main Zf -band source. We have re-measured the flux den- 
sities coming from the full extended nebula as seen on the 
UKIDSS images, after first removing all the interloping stars 
from the images and replace the pixels values using the adja- 




URL: http://www.leeds.ac.uk/RMS 



Fig. 4. Left: UKIDSS //-band (upper panel) and ZT-band (mid- 
dle panel) observations of the base of the W33A outflow scat- 
tering nebula. Images are logarithmically scaled, in units of 
MJy/sr and rotated by 35° so that the main lobe is pointing 
downwards, //-band contours are at 7.5%, 15%, 40% and 70% 
of peak flux; the Zf -band contours are at 1%, 2.5%, 10%, and 
40% of peak flux. Secondary stellar sources have been removed 
from the images. The lower panel shows the H — K colour im- 
ages (in magnitudes). Right: Model images produced with the 
RT scattering code, see Sect. |4.1.2l The model images are con- 
volved with 2-D Gaussian function with a FWHM of 1.6". All 
images cover an area of 25" on the side. 



cent regions. We find that in the Zf -band the nebula is 0.39 Jy 
(Zf ~ 8.0) and in the //-band 0.025 Jy (// ~ 1 1 .5). 

4. Modelling 

In this section we present the simultaneous modelling of 
the visibility and spectral energy distribution of W33A. 
Historically, the SEDs of MYSOs have been interpreted as hot 
stars embedded in a spherically symmetric dusty envelope (e.g. 
Giirtler & Henning 1986). In PapeiT we adopted the approach 
of modelling the observables with a 1-D spherical symmetric 
radiative transfer code. We concluded that models with shal- 
low density laws are preferred provided that a relatively low 
effective temperature for the central star is adopted. Here we 
construct a fiducial model of a hot star surrounded by a proto- 
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stellar envelope with outflow cavities. Our initial model aims 
to be as simple as possible and we begin by investigating the 
imprint on N-band visibilities by a protostellar envelope that in- 
cludes outflow cavities. The walls of the cavities are suggested 
to be the source of mid-IR emission seen on 1" scales in the 
nearly edge-on MYSO G35. 20-0.74 (De Buizer 2006) and in 
Ceph A HW2 (de Wit et al. 2009). 

4.1. Description of the radiative transfer code and 
input 

We adopt the 2D-axi-symmetric dust RT model by Whitney et 
al. (see for details Whitney et al. 2003a, 2003b). The model 
calculates radiation transfer (absorption, re-emission and scat- 
tering) through a dusty structure. The inner rim of the struc- 
ture is the dust sublimation radius which follows from a fit to 
models with different stellar temperatures (see Whitney et al. 
2004: R su bi/R* = (T su bi/T„y 2A ); space is empty within the in- 
ner rim except for the star itself. The structure can include three 
distinct geometrical elements: an accretion disk, a protostellar 
envelope, and low density polar outflow cavities. The protostel- 
lar envelope is described by the analytical TSC solution (Ulrich 
1976; Tereby, Shu & Cassen 1984) of a simultaneously rotating 
and collapsing spherical structure. This solution is governed by 
the infall parameter (Minfdi) and the centrifugal radius R c , i.e. 
the radius where rotational motion dominates infall. The enve- 
lope mass infall rate is a parameterization of the density dis- 
tribution and does not per se represent an actual determination 
of infalling material. The presence of an accretion disk is op- 
tional in the code. When present, the disk structure follows the 
one for a flared a-type disk prescription, and the released ac- 
cretion luminosity within the dust inner rim is added to the cen- 
tral luminosity source. The outflow cavity geometry can have 
two possible shapes, a streamline or polynomial geometry. The 
streamline is conical on large scales, the polynomial has a 180° 
opening angle at the stellar surface. Each of these geometrical 
elements are potential mid-IR emitters and could contribute to 
the visibilities. For most model runs the R sub i is comparatively 
small (tens of AU) and unresolved by the employed VLTI base- 
lines in the N-band. 

A model run is necessarily characterised by a large number 
of free parameters. A fraction of these can be constrained by 
the known properties of W33A. The spatial information at milli 
arcsecond angular resolution provided by MIDI allows the RT 
model to be seriously tested on its correctness in describing a 
MYSO. 

4.1 .1 . Fixed input parameters 

The central star parameters are chosen close to an 07.5 ZAMS 
according to the scale of Martins et al. (2005). The bolomet- 
ric luminosity is 10 5 L Q (Faundez et al. 2004) and it is lo- 
cated at the kinematical distance of 3.8 kpc (Bronfman et al. 
1996). The luminosity is based on IRAS flux density, and could 
be overestimated. We take the opportunity here to correct our 
70 pm Spitzer MIPS photometry presented in Paper I, which 
was not corrected for detector non-linearity. Taking this cor- 



rection into account we find a value of 1.9 10 3 Jy, quite similar 
to the IRAS measurement at 60 pm. This indicates that the lu- 
minosity based on IRAS is probably reasonable. We think it 
is unlikely that a significant fraction of the IR luminosity is 
due the secondary millimetre source noted in van der Tak et al. 
(2000) at 6" from the main source. The W33 A MIPS 24 pm and 
70 pm data show no sign of a secondary object, nor does a new 
24.5 pm image taken with Subaru/COMICS at an angular res- 
olution of 0.6" (in prep.). This suggests that the IR luminosity 
is dominated by a single object. 

A number of the important parameters that determine the 
properties of the TSC envelope are fixed for the following rea- 
sons. We choose a single type of dust that consists of an inter- 
stellar mixture of "warm silicates" (Ossenkopf et al. 1992), and 
amorphous carbon with a size distribution according to Mathis 
et al. (MRN; 1977). The type of silicates provides a better fit 
to the observed shape of the absoprtion feature than MRN DL 
grain-mixture (Draine & Lee 1984), in particular to the blue 
wing (see paper I, and see also Capps et al. 1978; Giirtler & 
Henning 1986). The centrifugal radius is set at 33 AU which 
is comparable to the dust sublimation radius. TSC envelopes 
make a transition from a pure infall, - 1 .5 radial density law to 
a shallower rotating, -0.5 law at the centrifugal radius. A TSC 
envelope with an centrifugal radius smaller than the sublima- 
tion radius is dominated by the infalling part of the envelope's 
structure solution, i.e. the envelope is quite spherically sym- 
metric apart from the outflow cavities. The shape of the outflow 
cavities is polynomial (see Whitney et al. 2003), which in the 
presence of a stellar wind could be a more appropriate geome- 
try than a conical streamline. The density in the outflow cavity 
is chosen to be typical of MYSO outflows, i.e. between 5 10 3 
and a few times 10 4 cm" 3 (see Beuther et al. 2002). The cavity 
material consists of gas and dust. We initially base the outer 
edge of the envelope on 350 pm observations by van der Tak 
et al. (2000). The radial intensity profile shows a downturn of 
submm emission at a radial distance of ~ 50", or 2 10 5 AU. We 
require a posteriori that the model densities at the outer edge 
are similar to that of the large scale molecular cloud material in 
which W33A is embedded, viz. 10 4 cm" 3 . 

4.1.2. Cavity opening angle 

One can exploit the observed shape of the near-IR scattering 
nebula to constrain the opening angle of the outflow. To this end 
we use the less computationally intensive scattering-only ver- 
sion of the dust RT code (Whitney & Hartmann 1992, 1993). 
The scattering code provides images at ten inclination angles p 
simultaneously, which one can readily compare to the observa- 
tions. We model the near-IR nebula under the assumption that 
the emission is free from direct thermal emission. We put par- 
ticular emphasis on matching the outflow opening angle with 
the observed one and constraining the inclination angle. Model 
input parameters are conform the fixed parameters discussed in 
the previous subsection, except for the dust properties. The de- 
tailed dust properties are not important for the goal to constrain 
the geometrical parameters of the W33A system from the near- 
IR nebula. For completeness, the scattering dust used here is 
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Table 2. Parameters for the preferred model discussed in Sect. 14. 21 This model consists of an envelope with outflow cavities; 
a disk is not present.. The inner radius was chosen equal to the dust sublimation radius. The model has a centrifugal radius of 
33 AU, and cavity opening angle of 29 = 20°. Most parameters are self-explanatory. M; nta n is the envelope's mass infall rate; i is 
the viewing angle with respect to the outflow axis; R env is the envelope's outer radius. 
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Fig. 5. Left: SED of the preferred model overplot W33A flux densities. Small filled circles trace the ISO SWS spectrum, small 
open circles are the MIDI spectrum, filled square is the Spitzer MIPS point, open square is the IRAS measurement at 60/mi. In 
the (sub)mm range, the integrated fluxes at 350^/m, 850 pm, and 1.2 mm are from van der Tak et al (2000), Di Francesco et al. 
(2008), and Faundez et al. (2004) respectively. Right: Zoom in onto the N band spectrum. Shown are the MIDI flux spectrum 
with errorbars and the model fluxes appropriate for the MIDI slit width of 0.6" 



coated with a small layer of water ice. The dust sizes follow 
a distribution that fits the extinction curve from Cardelli et al. 
(1989) for a prescription of Ry = 4.0. This distribution is not a 
single power law, and it is presented in Kim et al. (1994). For 
more details on the dust properties see Fig. 1 in Whitney et al. 
(2003a). 

Fig. 2] compares the H and /T-band images of a model run 
with the W33A scattering nebula. The images of the observed 
nebula have been rotated by 35°. The observed and model im- 
ages are in the same units (MJysr -1 ). Model images are con- 
volved with a Gaussian point spread function that simulates the 
one derived from the UKIDSS images and Gaussian noise de- 
rived from the observed images is added to the model ones. The 
outflow angle measured from the images directly would be 45°, 
however the presented model has an envelope geometry with a 
full cavity opening angle of less than half that, viz. 29 = 20°. In 
our fitting procedure, attention is devoted to the absence of an 
//-band counterpart of the central source. Whether or not a cen- 
tral source is detectable depends on the system inclination. The 
model that is presented in Fig.|4]has an inclination of 50°. An 
inclination less than 40° or more than 70° are excluded. In the 
former case the central star is visible in the //-band, whereas in 



the latter case the redshifted outflow lobe becomes prominent. 
These features are not present in the UKIDSS near-IR images. 

Although the shape on the sky can be reproduced reason- 
ably well, not a single model run reproduces the observed H—K 
colour of ~ 3.5. Photometry done on the model images shows 
that they produce colours that are much bluer instead. The 
model images of Fig. [4] match the colour of the observed im- 
ages, but additional foreground extinction of Ay w 25 is re- 
quired. The model H and Zf-band fluxes need to be scaled up 
accordingly. The central source visible in the /if -band consti- 
tutes only a minor contribution to the total observed flux den- 
sity at this wavelength. We note that MYSO models invariably 
have problems matching especially the short wavelength fluxes 
of MYSO scattering nebulae (e.g. Alvarez et al. 2004). The 
deviant model colours might be rectified by the inclusion of 
a clumpy structure within the outflow cavity (Indebetouw et 
al. 2006) and/or very small dust particles that are transiently 
heated to high temperatures. We find that the introduction of 
an outflow cavity with a smooth density distribution into an in- 
falling envelope cannot be the only physical element determin- 
ing the observed near-IR emission of MYSOs. Whichever the 
solution, a deeper investigation of this problem is beyond the 
purpose of the paper. Nonetheless, a comparison of near-IR im- 
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Fig. 6. Model (full line and filled circles) fits to the observed visibilities. Model parameters are listed in Table 2. Corresponding 
images are presented in Fig. [8] Note the discrepancy with the spherically symmetric model (dotted lines) used to fit configuration 
D in paper I. 



ages with scattering RT models has led to relatively stringent 
constraints on the outflow opening angle, and comparatively 
less stringent limits on the inclination. 

4.1.3. Methodology 

The prime model parameters which remain to be determined 
are the density as parametrized by the envelope mass infall 
rate, the presence of an accretion disk (and its properties), and 
the outer radius. They determine the peak of the SED, and 
depending on inclination angle, the depth of the silicate fea- 
ture, and the dominant emission regions on scales of 100 AU. 
The RT code is computationally demanding which excludes 
a grid search for the optimum value for the decisive model 
parameters. Building on a simple fiducial model, we try to 
find a model that fits the SED and MIDI visibilities simul- 
taneously. One run delivers the SEDs at 10 inclinations si- 



multaneously, evenly spaced in cos(/v). Model images are pro- 
duced for a single inclination at a time but for any number 
of desired wavelengths/filters. A total of 4 10 7 "photons" are 
used to produce images, for SED-only-runs the code uses four 
times less. We divided the A^-band up in 6 square filters of 
1 /im width, running from 7.5 //m to 12.5 jum. The resulting 
narrow-band model images are multiplied with a Gaussian 
model for the VLT-UT Airy disk. This represents the aper- 
ture field of view with a FWHM appropriate for the partic- 
ular wavelength. The resulting image is then Fourier trans- 
formed and the visibilities extracted for the projected interfer- 
ometer baseline and position angle. The model images have a 
pixel size of 0.01" which is less than half the angular reso- 
lution (0.024") at the longest baseline for the shortest wave- 
length. The SED is matched from ~8/zm to 1 mm. The input 
data are the MIDI flux spectrum, the slope of the SED around 
40yum determined from the ISO-SWS spectrum (see Paper I) 
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and the MIPS 70//m flux density. At the longer wavelengths, 
the integrated fluxes from the following observations are used: 
350 yum from with the CSO/SHARC instrument (van der Tak 
et al. 2000), 850/mi with JCMT/SCUBA (Di Francesco et al. 
2008), and 1 .2 mmm with SEST/SIMBA (Faundez et al. 2004). 
The far-IR and (sub)mm observations are compared to model 
SEDs representing all the emitted energy at each wavelength, 
whilst the MIDI flux spectrum is compared to the model SED 
with an aperture of 0.6" corresponding to the MIDI slit width. 
Matching of the SED and visibilities is performed by eye. 

4.2. Preferred model fit 

With the parameter values described in the previous section, a 
model is obtained that produces a reasonable match to both the 
SED and the visibilities. We have run a total of some 30 models 
in order to narrow down the parameter values varying the mass 
infallrate, inclination angle, and the outer radius. A satisfactory 
model is presented in Figs.[5]and[6] The parameter values of 
this particular model can be found in Table 2. The figures show 
that the model fits both the SED and visibilities reasonably. A 
discrepancy between model and observed SED is found in the 
near-IR, as expected (see Sect. 14.1.21 ). The model reproduces 
quite well the spectral trend of smaller visibilities in the wings 
of the silicate feature, especially clear on the blue wing. Fig. [6] 
also displays the visibilities from the spherically symmetric so- 
lution that matches the SED and visibilities of configuration 
D only (from Paper I). Each panel shows the spherically sym- 
metric model visibilities calculated for the appropriate baseline 
length. It is clear that they provide a considerably worse fit to 
the visibilities obtained on the additional configurations A, B 
and C, illustrating the need for multiple baselines. 

The actual model images of the preferred model are shown 
in Fig. [8] along with a slice through the spatial temperature and 
density distribution in Fig. [7] They make clear that the domi- 
nant emission regions in the iV-band at a few 100 AU scale are 
the envelope regions close to the surface of the outflow cavity. 
These envelope regions are warmed up due to irradiation by 
the star. We refer to these regions as "cavity walls" for brevity. 
The density distribution in Fig. [7] illustrates that the density of 
the paraboloidal cavity is a few order of magnitude less than 
that of the envelope. The cavity itself is hotter than the bulk 
of the envelope but contributes very little to the A^-band emis- 
sion due to the low density. Fig. [8] shows that the emission is 
relatively extended on the probed scales and produces conse- 
quently very low visibilities. An increase in optical depth, by 
either probing a higher inclination angle or increasing the mass 
infall rate, extinguishes the warmest, inner part of the irradi- 
ated walls. As a result visibilities generally decrease with an 
optical depth increase. However the extremely deep silicate ab- 
sorption feature of W33A requires a large optical depth. This is 
the trade off in order to obtain a decent fit. We find that for an 
inclination angle larger than 35° the mass infall rate needs to 
be smaller than 1 10~ 3 M yr _1 in order to see the warm cavity 
wall regions and avoid spatially resolving the A^-band emission 
out. The increase in optical depth due to the silicate absorp- 
tion has a similar effect. Fig.[8]makes clear that the warm inner 



regions of the cavity walls become progressively fainter with 
wavelength causing an increase of the typical emitting size and 
hence lower visibilities. 

The same effect is found when one increases the centrifu- 
gal radius. The preferred model with an increased R c to a value 
of 300 AU shows that the density in the equatorial region in- 
creases but the increase is confined to within a fraction of a 
degree of the equator. The density actually decreases away 
from the equator allowing a larger area of the cavity walls to 
be heated up. Visibilities drop therefore with increasing R c . 
Further fine tuning of the model parameters would require ad- 
ditional observations at preferentially shorter baselines. 

In summary, MIDI visibilities and SED of W33A can be 
matched with a simple geometry that consists of a protostellar 
envelope with a - 1 .5 radial density distribution that has polar 
outflow cavities with a 20 — 20° opening angle. We note that 
the modelling of the MIDI data is not very sensitive to the exact 
density powerlaw in the envelope within the framework of the 
applied RT model. The density powerlaw in this case is effec- 
tively governed by the centrifugal radius. Spherical models as 
applied to W33A in de Wit et al. (2007) are inadequate in mod- 
elling structure on 100 AU scales as Fig. [5] illustrates. On the 
other hand, resolved observations of MYSOs at ~ 25 jum at av- 
erage inclination angles (de Wit et al. 2009) and at 350 //m (van 
der Tak et al. 2000) provide insight into the density law starting 
at radii of 1000 AU and indicate a radial density law of -1.0. 
Matching images of W33A with model images at these respec- 
tive wavelengths will provide much stronger constraints on the 
centrifugal radius. For now we show in Fig.[9]the 350/im mor- 
phology of our model. We have convolved the model image 
with a 10" beam appopriate for comparison with the CSO ob- 
servations presented in van der Tak et al. (2000). The latter ob- 
servations show a structure more flattened than the model im- 
age, however the FWHM of the model image of 16" is brack- 
eted by the observed FWHM range of 14.4"-19.4", showing 
a good correspondance. The density at the outer radius of the 
model complies with expected densities of molecular clouds. 
Observations in the submm give a lower limit to the size, while 
in reality the outer parts of the envelope may become even 
steeper than -1.5 (e.g. Mueller et al 2002). The outer radius 
is therefore uncertain by a factor of 2. The best matching RT 
model presented does not need contributions to the emission 
from other components within the MYSO environment, instead 
iV-band emission on 100 AU scale is completely dominated by 
the warm cavity walls. 

5. Discussion 

In several studies the accretion based formation scenario for 
massive stars has been forwarded partially based on observa- 
tions with 1" angular resolution and/or SED studies of MYSOs 
(e.g. De Buizer et al. 2005). If MYSOs are to be considered 
young massive stars still in their accretion phase, and this ac- 
cretion occurs in a similar fashion as in low-mass stars, then 
circumstellar accretion disks ought to be present. The MYSO 
accretion disk is expected to be luminous (approaching that of 
the star itself), and estimated to have a size a few times 100 AU. 




Fig. 7. Logarithmically scaled slices through the spatial distribution of the density (left) and temperature (right) of the preferred 
model perpendicular to the outflow axis. The cavity has a very small, but non-zero density consisting of gas and dust (see text). 
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Fig. 8. Logarithmically scaled images (top row) and corresponding Fourier transforms (bottom row) of the preferred model 
(Table© at wavelengths 7.5 /mi to 12.5 /mi in steps of 1 jum. Image field of view is 0.6" (2280 AU) on a side, corresponding to 
the MIDI slit width. Pixel scale is 0.01". The model images have been multiplied with a wavelength dependent 2-D Gaussian 
function, simulating the VLT-UT airy disk. Warm cavity walls dominate the emission at 10/im. 



Our presented analysis prompts the question to what degree an 
accretion disk can contribute to the mid-IR emission. 



5.1. The contribution of an accretion disk 

The dust RT code allows the inclusion of a dust disk truncated 
at the 25 AU dust sublimation radius. As a first step we con- 
sider the effect on the model fit by adding such a dusty disk 
to the preferred model presented previously (Table 2). We em- 
phasize that the disk is directly irradiated by the star, as the 
space between disk and star is empty in the model. The mass of 
the dusty disk is initially chosen to be 1 % of the stellar mass, 
i.e. 0.25 M Q . The radial density law of the disk concurs with 
that expected of an a-type disk (i.e. a power law exponent of 



1.875, Frank et al. 1985). Due to the large inner truncation 
radius, the actual accretion luminosity generated by the dust 
disk is marginal and on the order of 1% of the stellar lumi- 
nosity (see eq. 6 in Whitney et al. 2003a). The formal accre- 
tion rate (see eq. 5 in Whitney et al. 2003a) for a = 0.01 is 
then 2.8 10~ 6 (a/0.01) M yr _1 . We note that the accretion lu- 
minosity inside the dust sublimation radius, where no actual 
disk structure is present in the RT code, is emitted with the 
stellar spectrum (see Whitney et al. 2003a). The disk has a con- 
stant opening angle, i.e. the scaleheight increases linearly with 
radius. The model disk outer radius is arbitrary because it does 
not emit in the A^-band, and is better constrained by millimetre 
data; for now, it is set to 500 AU. 
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Fig. 9. Morphology at 350 /mi predicted by the preferred 
model. Contours are at 15%, 30%, 50%, 70% and 90% of the 
peak value. Model image has been convolved with a 10" beam. 
See van der Tak et al. (2000) for the observed 350 //m image. 



The resulting model visibilities of a protostellar envelope 
plus dust disk are shown in Fig. [10] It is clear that such a model 
violates the observed MIDI visibilities shortward of the sili- 
cate absorption feature. At these wavelengths the dust disk con- 
tributes ~ 50% to the total flux density, and one can deduce that 
the disk has a visibility of approximately 70%. In order for this 
model to be consistent with the observations, we compute that 
the fractional contribution to the total 8 /mi flux by the dust 
disk should be at most 1%, if the disk emission at this wave- 
length is fully unresolved. If it is marginally resolved (as it is) 
the 1 % upper limit becomes smaller. By running models with 
increasingly smaller disk masses, we find that this flux frac- 
tion is reached for a total mass of the dusty disk of 0.01 M Q , or 
equivalently 1.1 lO-^a/O.O^Moyr 1 . The SED of W33A is 
still fitted well by the dust disk+envelope model. We thus find 
that the disk mass is very small compared to the envelope mass, 
and for reasonable values of a that the dust disk accretion rate 
is more than three orders of magnitude smaller than the enve- 
lope mass infall rate. These results on the disk masses are in 
accordance with other studies of MYSOs (e.g. Murakawa et al. 
2008) and also with some evolved early-type Herbig Be stars 
(Alonso-Albi et al. 2009). Would the disk outer radius be much 
smaller than the assumed 500 AU, and similar to the small outer 
radii found for the Herbig Be stars, then the disk mass should 
be scaled down accordingly; a 50 AU outer radius results in 
a disk mass of 5 10~ 4 M G . A small disk mass may imply that 
such a dust disk is not present, or that it is simply extremely 
faint in the mid-IR due to the absence of gas and dust opacities 
in the model disk, which would lead to strong shielding from 
the stellar radiation field. 

Most of the accretion luminosity however would be re- 
leased within the 25 AU dust sublimation radius by the hot 
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Fig. 10. Model predictions for baseline configuration A. 
Represented are the preferred model (full line), with the addi- 
tion of a dust disk (M - 0.25 M , dashed line), and an optically 
thick disk interior to the dust sublimation radius (dotted line). 



gaseous extension of the dust disk. Due to self-shielding, such a 
disk probably contains dust at smaller radii than the 25 AU sub- 
limation radius determined from direct irradiation. Inclusion 
of a self-luminous accretion disk within the dust sublimation 
radius is required to fit the SED and the near-IR interfero- 
metric data of the Herbig B6e star MWC 147 (Kraus et al. 
2008). Accretion luminosity generated by an accretion rate of 
7 10~ 6 M o yr _I for this particular object is however not essen- 
tial to model the star's mid-IR interferometry, which is already 
reasonably well reproduced by a passive disk. We explore to 
what extent the mid-IR interferometry of W33A is affected 
by the inclusion of an optically thick accretion disk. We ap- 
proximate the reslulting fluxes and visibilities with an optically 
thick, geometrically thin a disk model (see Malbet & Bertout 
1995; Malbet et al. 2007) and add these to the preferred en- 
velope model in Table 2. The disk is heated by stationary ac- 
cretion only and for an accretion rate equalling the mass in- 
fall rate (7.5 10 _4 M Q yr _1 ) the generated accretion luminosity 
is ~ 7 10 4 L for our adopted central star. The spatial properties 
of the optically thick disk emission is such that 99% of the N- 
band is emitted within 50 AU. Most of the 8 /mi flux is emitted 
at 20 AU. The disk thus remains largely unresolved on the em- 
ployed VLTI baselines with MIDI. Disk accretion at a slower 
rate than the one assumed are cooler and thus more compact. 
We choose the inner radius of the disk to be at 1 R,. The disk 
temperature falls below 1500 K at 2.5 AU, beyond which the 
disk opacity is given by both gas and dust. 

As in the previous paragraph, if the preferred model visi- 
bilities are to be affected by the optically thick disk, the ratio of 
disk emission to the emerging envelope emission should be at 
least 1 % at 8 /mi. We calculate the resulting disk spectrum for 
the said disk accretion rate and apply an extinction of Ay=230 
due to the overlying envelope, as found in the line-of-sight to 
W33A from the preferred model. We show the resulting visibil- 
ities in Fig. [10] The figure shows that an a disk model embed- 
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ded in a protostellar envelope is compatible with the observed 
MIDI visibilities. The model SED does not differ from the en- 
velope only SED, provided that the total luminosity of each 
model is chosen the same. The a disk contributes less than 
1% to the total 8//mflux, despite its high luminosity. The en- 
velope of W33A is simply too opaque for any disk emission to 
be directly observed. MIDI observations of MYSOs with less 
massive envelopes or at smaller inclinations however may well 
disclose the presence of an accretion disk in such systems. The 
expected visibilities would be much larger than the ones ob- 
served for W33A (depending on distance and accretion rate), 
considering the compactness and brightness of the 8 fim emis- 
sion in accretion disks. 

The popular YSO SED model grid developed by Robitaille 
et al. (2006J1 allows an estimate of the extent to which MYSO 
SEDs are dominated by dust disk emission. This YSO SED 
grid makes use of the same 2-D axi-symmetric RT calculations 
as we use in this paper (Whitney et al. 2003), and the disks 
in the SED model grid are therefore dust disks. We extract all 
(135) SEDs corresponding to ZAMS hot stars with a luminos- 
ity of 10 3 < L/Lq < 10 6 and an envelope mass infall rate of at 
least 10~ 4 M yr~' . Inspection of the fractional disk flux contri- 
bution at 8 fim (just short of the extended wings of the silicate 
absorption feature) shows that all except three of the extracted 
models have a fractional contribution larger than 15%, and half 
of the models have a contribution larger than 70%. The geom- 
etry of the accretion disk determines the total visibility at these 
wavelengths. The SED model grid thus contains a preponder- 
ance of MYSO models with relatively high visibilities, simi- 
lar to the envelope plus dust disk model presented in FigfTOl 
At present, the MIDI observations of MYSOs demonstrate the 
opposite trend however. Visibilities presented for W33A are in 
line with the low /V-band visibilities observed for other MYSOs 
(see Linz et al. 2008; Vehoff et al. 2008; de Wit et al. in prep.). 
The latter could equally well be explained in the context of cav- 
ity wall emission as proposed here for W33A. If verified, then 
MIDI has not (yet) revealed the presence of a MYSO accretion 
disk on size scales of ~ 100 AU, but has certainly the potential 
to do so. 

5.2. SED model grid fits 

To what extent are models obtained from SED fits consistent 
with the spatial information on milli arcsecond scales? To ad- 
dress this question, we apply the web-based SED fitter proce- 
dure developed by Robitaille et al. (2007) to the W33A SED 
data shown in Fig. [12] The procedure returns a ranking on the 
basis of a x 1 minimization method, and we discuss the grid 
models following this ranking. 

The procedure delivers as best-fit a model consisting of a 
central hot star and a dust disk (model grid number 3002520, 
see Table|3). We see from Fig.[T2]that the grid model fit to the 
SED is acceptable, as expected. This particular grid model has 
a system inclination of 50°, consistent with the one of the re- 
flection nebula. The best-fit model consists of a protostellar en- 
velope with a cavity opening angle relatively large compared to 

2 URL: http://caravan.astro.wisc.edu/protostars/index.php 



the one derived from the scattering nebula. Moreover a fit of the 
model to the observed SED requires a distance of 1 kpc, which 
is a factor ~ 4 too small. The discrepancy between the modelled 
and observed SED is mainly at the depth of the silicate feature. 
This depth is better accounted by dust partially consisting of 
"warm silicates" (see Ossenkopf et al. 1992) as discussed pre- 
viously; the appropriate set of optical constants is not used in 
the SED model grid. Broken down in its various contributing 
components, we find that the grid model SED is completely 
dominated by disk emission at wavelengths shortward of the 
silicate absorption feature. The emission accounts for the rel- 
atively high flux levels observed shortward of the silicate fea- 
ture. This is a common property of MYSO SEDs. Spherical 
and indeed 2D-axisymmetric RT codes for the protostellar en- 
velope are unable to reproduce these high flux levels (see e.g. 
de Wit et al. 2009). Contrary to the good SED fit, Fig.[T3lshows 
that the corresponding visibilities are far too large. The visibil- 
ities are obtained by producing images with the dust RT code 
using the same parameter file that produced the SED. The grid 
model visibilities are inconsistent simply because of the domi- 
nance of the (compact) disk emission. Clearly, a good model fit 
to the SED does not guarantee a proper description of the ob- 
ject in consideration; the best-fit grid model violates the spatial 
information delivered by MIDI. 



If we use the known properties of W33A and limit the 
model distance to be between 3.0 and 4.6 kpc with a maxi- 
mum foreground extinction of A v = 50, then we find that the 
returned fitting models have a strong preference for small incli- 
nations. In the best 10 models, 7 have an inclination of 18.19 
degrees (the smallest present within the grid), also the best fit 
grid model (number 3004217). In addition, the visibilities fit 
the MIDI observations reasonably well (Fig.fTSl. The model 
is characterized by a central cool star with a large stellar ra- 
dius. Such a description of the MYSO M8E-IR is preferred 
in Linz et al. (2009). In the case of W33A, the small inclina- 
tion angles are inconsistent with the allowed range we derived 
from the scattering nebula. The predicted /f-band morphology 
of this model is shown in Fig.QT| and is strongly dominated by 
the central source which is in contradiction with the observa- 
tions of W33A (Fig. [4]). The ranked fourth model however has 
a more reasonable inclination angle of 50° (number 3009737). 
The latter model consists of a hot star and an accretion disk and 
encounters the previously discussed problems in fitting the vis- 
ibilities (Fig.fTJb. If we would apply in selecting a model from 
the grid all the W33A fiducial parameters (40° < i < 70°, cav- 
ity opening angle 10 < 20 < 30) for a hot star (T efi > 20 000 K) 
without a dust disk then the model complying with these con- 
straints is ranked 245th (model number 3003047). 



It is clear that the SED model grid can deliver models that 
fit the SED, but all the SED fitting models presented in this sec- 
tion violate the MIDI visibilities. Regarding W33A, the SED 
grid does not provide a fit to both SED and visibilities simply 
because the appropriate models are not present in the grid (see 
also the discussion in Robitaille 2008; Linz et al. 2009). 
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Table 3. Models from the YSO SED model grid (Robitaille et al. 2006), obtained using the web-based fit procedure. The quoted 
distance is required to fit the SED of W33A; is the outflow cavity opening angle; A v is the applied foreground extinction. 
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Fig. 11. Logarithmically scaled /T-band image for grid model 
3004217 to be compared with the observations in Fig. [4] The 
contours are at 1%, 40%, 50% and 65% of the peak flux. 
The image is convolved with a 2-D Gaussian function with a 
FWHMof 1.6". 



6. Conclusions 

Observations with MIDI at the VLTI of the MYSO W33A indi- 
cate that A^-band emission on scales of 100 AU is due to warm 
parts of the envelope close to the outflow cavity, i.e. the "cavity 
walls". The overlying envelope, when described by a TSC type 
geometry is characterised by an equivalent mass infall rate of 
order 10~ 3 M Q yr _! and is quite massive. We find that the pres- 
ence of a dust disk located beyond the dust sublimation radius 
violates the observed visibilities. In order to make the model 
consistent with the observations the dust disk is required to be 
marginal in mass. In contrast, the addition of an optically thick 
accretion disk interior to the dust sublimation radius has little 
effect on the visibilities, even when the disk dominates the lu- 
minosity. The disk emission is highly extincted by the massive 
overlying protostellar envelope. MYSOs viewed at smaller in- 
clination and/or with less massive protostellar envelopes have 
the potential to reveal such accretion disks. 




Fig. 12. Predicted SEDs of the discussed models obtained from 
the SED web fit procedure listed in Table[3] 
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